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Abstract 

We investigate the question of optimal input ensembles for memory channels 
and construct a rather large class of Pauli channels with correlated noise which 
can be studied analytically with regard to the entanglement of their optimal input 
ensembles. In a more detailed study of a subclass of these channels, the complete 
phase diagram of the two-qubit channel, which shows three distinct phases is 
obtained. While increasing the correlation generally changes the optimal state 
from separable to maximally entangled states, this is done via an intermediate 
region where both separable and maximally entangled states are optimal. A more 
concrete model, based on random rotations of the error operators which mimic 
the behavior of this subclass of channels is also presented. 
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1 Introduction 



A basic question in quantum information theory [U |2l El [U [5] is whether the use of 
entangled states for encoding classical information can increase the rate of information 
transmission though a channel or not. A proper calculation of the so called Holevo 
capacity [3j of a channel, representing by a Completely Positive Trace preserving 
(CPT) map $, requires the optimization of Holevo information over ensemble of input 
states when we encode information into arbitrary long strings of quantum states (more 
precisely states in the tensor product of the Hilbert space of one state) and carrying 
out the limiting procedure C := lim„_»oo C„,, where 

Cn ■■= -SupeXnie) (1) 

n 

is the capacity of the channel, when we send strings of n quantum states into the 
channel. Here e := {pi,Pi} is the ensemble of input states, 

Xn{e) := SiY^pMPr)) - Y.P^Si'^(Pi)) (2) 

i i 

is the Holevo information of the ensemble and S{p) = —tr{p\ogp) is the von Neu- 
mann entropy of a state p. To find the capacity of a given channel we should find 
the ensemble which maximizes this quantity and we call it the optimal ensemble of 
input states. Then the properties of this ensemble can be studied and ask whether 
this ensemble includes entangled states or not. 

The importance of this question stems from the fact that entanglement is a vital quan- 
tum mechanical resource in many tasks in quantum information processing, however 
it is a held belief that entanglement is so fragile in the presence of noise. So it would 
be interesting if one can show by using this property higher rate of data transmission 
through noisy channels is achievable. The difficulty in answering this question is not 
only because of the optimization of Holevo information over multi-parameter space 
but also due to the fact that no concrete classification of multi-particle entangled 
states exist. 

A much simpler problem is to calculate C2, rather than C„, which is equivalent 
to refresh the channel after each two uses, and see whether entangled states can en- 
hance Holevo information or not. This problem has been tackled by many authors 
P El El [ini [III [la [la [111 [ig [16] and for a kind of correlated channel it has been 
shown that entanglement can enhance the Holevo information, if the correlation is 
above a certain threshold value. This correlated channel was first introduced in [9] as 
follows 

3 

^{p) = X! ^ CTjpai aj, (3) 

i,j=Q 

where ai, i = 0, 1, 2, 3 represent Pauli operators /, crx,(Ty , cr^ and Pij, the probabilities 
of errors cxj aj are correlated in a special way, namely 

Pij = {I- p)piPj + p6ijPj. (4) 
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The parameter // G [0,1], called memory factor signifies the amount of correlation 
in the noise of the channel. For /i = 0, the errors on the two consecutive qubits are 
completely un-correlated, while for /i = 1, the two errors are exactly the same. The 
idea behind this, is that when the channel relaxation time is much less than the time 
interval between the passage of the two qubits, the errors will be the same. What was 
observed in p] and in subsequent works [HI [HI [16] was that when fj, passes a critical 
value fic, the optimal input state jumps suddenly from product states to maximally 
entangled states. 

It is important to note that although the sharp transition of optimal ensemble of 
input states from product to maximally entangled states is interesting and cause non- 
analytical behavior of the channel, it is not obvious that the same happens for all 
kinds of correlated channels, the correlation ^ is one out of many different forms 
of correlations that one can envisage for a correlated noisy channel. Apart from the 
mathematical possibility of defining many other forms of correlations, one can also 
argue on physical grounds, in favor of other forms of correlation. A fully correlated 
channel which exerts a noise operator on the first qubit, need not exerts the same 
error operator on the second qubit, as in ([4]). In fact it is natural to expects that on 
exerting the first error, the state of the environment will change and depending on 
this new state, it will exert new errors with conditional probabilities on the second 
qubit even if the time interval between two consecutive uses of the channel be small. 

In this paper we try to shed light on these issues and to provide a basis for study- 
ing more examples and extend the study of correlated noisy channels by considering 
more general forms of correlations. Our study will be along the line of references 
[9l \T2[ [T^ 116] , that is we do not consider a specific model of environment, rather we 
take the abstract definition of the channel as a CPT linear map, defined by its Kraus 
decomposition [17| . We start in section 2 with the most general form of the action of a 
correlated Pauli channel on two qubits, where by two plausible requirements, we will 
restrict the parameters so that the members of this class can be studied by analytical 
means. Then we focus on one particular subclass and study in detail the optimal 
ensemble which maximizes its Holevo information. In this same section we propose a 
general definition of the correlation parameter, and show that when the correlation 
of the noise in this channel increases, the optimal ensemble changes from a product 
ensemble to a maximally entangled one. The phase diagram of the model also shows 
other interesting transitions not already observed in other works [9l [T2l [T^ [T6] . In fact 
the phase diagram contains three distinct phases, a product one and two maximally 
entangled ones, which differ with each other by a relative phase. We also observe that 
when we consider the correlation parameter and move through this phase diagram, 
the optimal ensemble changes from product to maximally entangled, however this 
transition is mediated by a region of correlation, where both maximally entangled 
and separable ensembles are optimal. Finally we construct a rather concrete model 
for this particular type of correlation. The paper concludes with a discussion. 
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2 The correlated action of a Pauli channel on two qubits 



The general action of a Pauli channel on two qubits is defined by the following CPT 
map 

^(p) =Y1 ^ij^i ^ ^jP^i ^ ^i' (^) 

where Pij is the probability of the errors aj and o"j on the first and second qubits 
entering the channel respectively. The number of independent error probabilities are 
15 due to normalization J2i,j Pij = 1- 

In order to see how much correlated the noise on the two consecutive qubits are, 
we form the marginal error probabilities p^^^ and p^"^^ and evaluate the following dis- 
tance between the two probability distributions, pij and pf''pf\ denoted as C 

where p^P := YlijPij ^'^'^ pT^ = J2j Pji- For uncorrelated noise we will have C = 
and for a fully correlated noise C will attain its maximum value. For the probability 
distribution C turns out to be C = pJ2i=oPi{^ ~Pi)i hence for fixed error proba- 
bilities, the correlations indeed increase with fi. 

There are two difficulties in studying such a channel for obtaining its optimal in- 
put ensemble and understanding how it depends on the correlation of the noise. The 
first one is that the probability distribution Pij can be correlated in many different 
ways and picking out a single parameter and designate it as the memory of the chan- 
nel, is just one possibility. The second problem is that the manifold of input states 
has itself 6 real parameters which means that the optimization task should be done 
over a 6 parameter space and thus analytical treatment of such a channel is almost 
impossible. To overcome these problems, we impose the following symmetry on the 
channel, 

^(c^3 ® (g) 0-3) = <!>{p). (7) 



This is the symmetry which has been considered in |T2] for making the model amenable 
to analytical treatment. Demanding this symmetry reduces the number of parameters 
in Pij to 7 parameters: 
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with normalization relation between the parameters. However it is more plausible to 
assume that the marginal error probabilities on the first and the second qubits be 
equal, that is 
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Here we are assuming that the errors on a sequence of qubits should be the same, 
regardless of how we enumerate the qubits of the sequence. What really matters is 
that any two consecutive uses of the channel are correlated. 



This assumptions in addition to the constrained composed by the symmetry in ([7] 
reduce the number of parameters from to 6 and the final form of the matrix of prob- 
abilities P with elements can be parameterized as follows: 
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where p + q + r + r] + s= ^. 

The advantage of demanding the symmetry in ([7|) in not only in reducing the param- 
eters of P but also in reducing the parameters of the general input states. 



Following [12j, we form the optimal ensemble by finding a state p* which minimizes 
the output entropy and hence minimizing the second term in the right hand side 
of ([2]). The input ensemble is then formed as a uniform distribution of the states 
£ = {pij := ((Tj aj)p*{ai ® The reason is that the first term of the Holevo 

quantity is maximized by this choise, i.e. 

= S{\l) = 2, (11) 

where in the first line we have used the covariance property of the channel: 

^{ai ajpai O aj) = {ai (g) aj)^{p){ai aj). (12) 

and in the second line the Shur's first lemma for irreducible representations of the 
Pauli group. 



Therefore finding the optimum ensemble of input states reduces to finding a sin- 
gle input state which minimizes the output entropy and we call it optimum input 
state. Regarding the convexity of entropy we deduce that we should search for the 
optimal input state among pure states which in general has 6 parameters. However, 
we are able to restrict the form of the input states to the simple states which are 
invariant under the above symmetry in ([7]) 

IV') =cos0|OO) +sin0e*'^|ll), (13) 

It is easy to find the output state of this channel with the above form of error probabil- 
ities. A straightforward calculation shows that the output state, in the computational 



4 



basis {|00), |01), |10), |11)} will be 
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Due to the block diagonal structure of this matrix, its eigenvalues can be calculated 
in closed form and hence a complete specification of the optimal input states can be 
made for this general class of correlated Pauli channels. 

2.1 Detailed study of a subclass and its phase diagram 

As an interesting and simple example, we consider a channel where only the param- 
eters p, q and r are non- vanishing, that is matrix of probabilities has the following 
form 

/ p \ 
g r 
r g 
\ p / 

where due to normalization p + q + r = ^ . The action of the Pauli channel on two 
consecutive qubits will then be: 



P 



(16) 



pp + pa^ ® (Tspcrs <8) 0-3 + qai (g) aipai ® (Ti + qa2 (8) o"2/9(72 (8) (T2 



+ rUi (g) (72/9fTi (g) (T2 + rcr2 (8) (71/9(72 (8) (Tl. 

The correlation parameter of this channel is found to be 

Using (fT5]) , the eigenvalues of the output density matrix are 

1 
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where 

Y := q{r — p) cos^ (j) + r{q — p) sin^ (f). 

For minimization of output entropy which is S{^{p)) = — AilogAi — A2logA2, we 
should maximize the difference between Ai and A2. This is achieved by minimizing 
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Figure 1: Color Online. The phase diagram of model (14). In each phase the minimum 
output entropy state (un- normalized) is specified. The triple point is {q,p) = (g, g)- 

Ysin^26. Therefore if y < 0, we should choose 9 = j and if y > 0, we should 
choose 9 = 0. Therefore the line Y = determines the boundary of the maximally 
entangled (6* = f ) and the product {9 = 0) optimal states. This line is determined 
by the equations q{r — p) = and r(q — p) =0. In the {q, r) plane these two lines are 
specified by the equations 2q + r = ^ and 2r + q = ^. 

When 9 = 0, the value of Y is immaterial, however when 9 = then the value of 
Y should be minimized again. In the maximally entangled region, where 9 = j, we 
should still minimize Y. For q{r—p) > r{q—p), Y takes its minimum at (p = 0, and for 
q{r — p) < r{q — p) it will take its minimum at (j) = ^. Thus the line q{r — p) = r{q — p) 
(q = r) will separate two types of optimal maximally entangled states from each other. 
The phase diagram is shown in figure ([T]). 

It is seen from the phase diagram ([2]) that as the correlation parameter increases, 
the optimal input ensemble changes from product to maximally entangled. There are 
however two remarkable features in this diagram not encountered in previous studies. 

First we see that depending on the values of the channel parameters, two different 
types of maximally entangled states, namely :^(|00) + |11)) and :^(|00) +i|ll)) are 
optimal. Although these two types of states, are transformed to each other by a local 

operator I ^ J ^ ) ' ^^'"^ channel is not covariant under this local operator, 

they should be considered different as far as optimality of the encoding is concerned, 
although they are equivalent as far as their entanglement properties are concerned |20j . 

Second, when we draw the contours of constant correlations in this phase diagram. 
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Figure 2: Color Online. The phase diagram of model (14). In each phase the min- 
imum output entropy state (un-normalized) is specified. The contours of constant 
correlations are also shown. As the correlation C increases, the optimal ensemble 
changes from separable to maximally entangled one. For ~ 0.43 < C <w 0.5, both 
types of states are optimal, the grey region shown in the inset. 

we observed that there is a region of correlation, for which both separable and maxi- 
mally entangled states are optimal, depending on the values of the parameters q and 
r, figure ([2]). We say that the two phases coexist, a property which is reminiscent 
of first-order transitions. We should stress that if one draws the phase diagram of 
the model in [12], not in terms of the parameter /i, but in terms of the correlation 
parameter C, one will again see such coexistence region. Therefore and specially 
with regard to recent studies in relating transitions in channel capacity to the critical 
transitions in their environment [181 119] , it is an interesting issue to see if transitions 
in the channel should be characterized as first or second order. 



2.2 A concrete and intuitive model of correlation 

In this section we construct a particular model for correlated noise in Pauli channels, 
which will reproduce the above example of correlation in a natural way. Consider a 
noisy Pauli channel acting on a qubit, defined as 
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i=0 
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where pi is the probabihty of error cjj (ctq = /) and J2i=oPi = 1- When the first qubit 
passes through the channel, and an error operator cjj acts on it, we assume that the 
state of the channel changes randomly and therefore on the second qubit, it exerts 
not the same error or a fixed error for that matter, but a random rotation of the fjj 
operator, in the form 

■■= U^fiTiU^g, (22) 

where U^^q is a random rotation around the axis n with angle 9. Thus a has the effect 
of the first error operator and also the random change in the environment. This ran- 
domness in contrast to a deterministic change in the environment state is physically 
plausible in view of the macroscopic nature of the environment. 

Therefore the action of the channel on two consecutive qubits may be written as 
follows 

3 

^{p) = Pop + X! Pi^i p ai(^ai. (23) 

i=l 

Since the rotations are random, the complete definition of the channel will be given by 
integrating over the above action with a suitable probability distribution over random 
rotations. Thus the final definition will be 

^„{p) = J diideP{n,0)^{p). (24) 

Clearly one can add more parameters to the above model, for example by taking 
different rotations to be along different axes. For simplicity, let us restrict ourself 
to a simple example in which the direction of all rotations are fixed in the z- axis 
and only the angle of rotation is random. Also to ensure the symmetry ^ we take 
Po = P3, and pi = P2, where po + Pi = ^- We take the probability distribution to be 
a Gaussian with mean value 9o and variance a. Hence the channel will be defined by 

^„{p) = / e-^^$(p), (25) 

where in this case di = e~2^'^^Uje2^'^^. One can say that parameter a is related to the 
memory of the channel. When o" = 0, the channel has full memory and it will exert a 
definite error operator (exactly the same error in the case = 0) on the second qubit 
depending on the operator which it has exerted on the first. However for a non-zero 
small value of o", the channel exerts errors on the second qubit which are close to 
the errors on the first qubit. As a increases further the memory is lost further and 
the channel will exert errors from a larger neighborhood of the errors on the first qubit. 

A remark is in order about the Gaussian distribution. The rotation operators are 
periodic which restrict the range of integration of 6 to [0,27r]. However this makes 
the subsequent formulas unduly cumbersome without adding much to the physics. 
Instead we can assume the variance a to be sufficiently less than 27r so that we can 
safely extend the range of integration of to (— oo, oo) and use the simple results of 



8 




0.2 0.25 0.3 0.35 0.4 0.45 0.5 



Figure 3: Color Online. The phase diagram of model (24). In each phase the minimum 
output entropy state is specified. 



Gaussian integration. 

After rearranging and doing the integrals, one finds that this channel has the form 
(jl7p with the parameters as given below 

P = Po, (26) 
1 + 

Q = Pii—-^ ), (27) 

r = Pii ^~l ), (28) 

The two independent parameters of this channel can be taken to be pi and a. In 
terms of these parameters the phase diagram is shown in figure Since we have 
always r < q, the region with "^(|00) + i|ll)) optimal state is not covered in this 
new phase diagram. The line which separates the product phase from the maximally 
entangled phase is now given by ()29p . Inserting the values of q and r from (j27p and 
(j28p in the relation 2r + q = ^ and simplifying we obtain 

e-2-' = 3-— . (29) 
Pi 

The phase diagram in figure© is re-drawn in terms of the new parameters in figure 
©• 

It is seen that depending on the value of pi , the optimal ensemble changes from 
separable to maximally entangled phase, when the memory passes a certain threshold 
(note that here a lower value of a means a larger value of memory). Also there are 
values of pi, where the optimal ensemble is always a maximally entangled one, no 
matter how weak the memory is. This is related to the fact that for no value of the 
parameter a, this channel is a product channel. 
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As stressed in [16] . the effect of memory on the type of optimal ensemble input 
and ultimately on the capacity of a channel, can be decided only when one considers 
the action of the channel on arbitrary long sequences of (entangled) qubits. With 
the type of memory introduced in [9], such extension is very difficult to pursue ana- 
lytically and one can consider very limited class of states. However with the type of 
memory introduced above, we think that such an extension is indeed more tractable 
by analytical means. For example one can consider the following type of correlated 
action of a Pauli channel on strings of n qubits: 

^p) = I dUidU2 ■ ■ ■ dUn-iJ2pM"\p) (30) 

i=0 

where 

^^r\p) = K^""'^ ® • • • ® <^M^t'^ ® • • • -f^ ® -^]\ (31) 

(k) t 

and (Tj = UkC^iUl in which Uk is a random unitary operator. Taking the random uni- 
taries from a distribution, one may relate the memory of the channel in a qualitative 
sense to the properties of the distribution as we did above. 

3 Discussion 

We have considered a large class of two-qubit correlated Pauli channels for which 
the entropy of the output state can be determined in closed analytical form. These 
channels are covariant under the action of the two-qubit Pauli Group and have the 
symmetry ^{p) = ^{(Tz'SKTzpcTz^crz)- For a subclass of these channels we have explic- 
itly determined the optimal input ensemble which maximizes the Holevo capacity and 
have determined the exact phase diagram, showing different regions where separable 
or maximally entangled states are optimal. There are two new features of this phase 
diagram. First, there are three phases separable and two different types of maximally 
entangled states which are optimal. Second if we use a precise definition of corre- 
lation, then as the correlation increases the type of optimal state generally changes 
from separable to maximally entangled states, however this is done through an in- 
termediate region where both separable and maximally entangled states are optimal. 
Put differently, if we think of correlation as a control parameter, the transition in our 
model is reminiscent of first order transition where there are regions of coexistence 
of the two phases. An interesting question is whether such transition occurs when 
the channel acts upon a string of qubits, not only two consecutive qubits [16]. This 
is a question which should be addressed if we are to judge definitely wether or not 
encoding of classical information enhances the classical capacity of quantum channels. 
Unfortunately settling this question seems to be very difficult. 

Finally one would like to see if the transition in the quantum capacity of quantum 
channels can be related to the well-studied subject of critical phenomena. In this 
direction a formalism has been developed in |18[ [T9] where models of correlated noise 
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are constructed by taking a many body system as the environment. The qubits which 
pass through the channel interact with this many body environment and the corre- 
lations in the many body state gives rise to correlations in the noise in the channel. 
The question of whether the transitions in channel capacity are of first order or not 
certainly have relevance for such studies. 
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